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ABSTRACT. In this review article, we present a unified approach to solving dis- 
crete, integrable, possibly non-commutative, dynamical systems, including the Q- 
and T-systems based on A r . The initial data of the systems are seen as cluster vari- 
ables in a suitable cluster algebra, and may evolve by local mutations. We show that 
the solutions are always expressed as Laurent polynomials of the initial data with 
non-negative integer coefficients. This is done by reformulating the mutations of 
initial data as local rearrangements of continued fractions generating some particu- 
lar solutions, that preserve manifest positivity. We also show how these techniques 
apply as well to non-commutative settings. 
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1. Introduction 

1.1. The recursions. In these notes, we mainly deal with the so-called Q- and T- 
systems based on the Lie algebra A r . These systems and their generalization to other 
Lie algebras were introduced [20] [22] in the combinatorial study of the complete- 
ness of the Bethe Ansatz states for the diagonalization of generalized Heisenberg 
quantum spin chains. Apart from their original representation-theoretic interpreta- 
tion as recursion relations for Kirillov-Reshetikhin characters (Q-systems 11201 ) or 
^-characters (T-systems |[T4l IT251 ), these equations have reappeared since in various 
combinatorial contexts lfl8l ED (26). We use here harmlessly renormalized vari- 
ables that allow, by introducing sign flips, to connect the systems to cluster algebras 
with trivial coefficients, and by a slight abuse of language, we will still call the cor- 
responding equations Q- and T-systems. 

Let I r = {1, 2, r}, we define the A r <2-system: 

(1-1) Ra,k+\Ra,k-\ = R 2 ai k+Ra+\,kRa-\,k {(XeI r ;keZ) 

(1.2) R 0jk = R r+hk = l (*€Z) 

and the A r T-system: 

(1-3) T a j£+\T a ,j,k-l = Ta^^i^Taj-i^ + Ta+ij^Ta-ij^ (dt G I r ',j,k € Z) 
(1-4) T j, k = T r+1Jjk =l 
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This involves quantities T a j & with fixed parity of a + j + k (which we take = mod 
2). 

The latter system is readily seen to reduce to the former, upon "forgetting" the 
index j. One way to realize this is to view the solutions of the ^-system as those 
of the T-system that are 2-periodic in the index j, with the correspondence R a ,k = 

Ta,a+kmod2,k- 

The T-system may also be viewed as a non-commutative generalization of 
the <2-system, in which R a n — > Rq. „ is no longer scalar, but an invertible operator 
acting on a Hilbert space with basis {\j)}jez- In this formulation, T a j^ are the 
matrix elements of the operator R a ^ in the basis. For a = 0, 1, r+ 1 we define the 
actions: 

(1.5) ^a,k\j + k+a)=T ajjjk \j-k-a) (j,k£Z) 

Defining also the shift operator d acting as d| j) = \j — 1) and its formal inverse d" 1 , 
the T-system is readily seen to be equivalent to the non-commutative A r Q-system: 

(1.6) R^+iR^Rca-i = Ra,*+R« + i,*R-iRa-i,* (aeI r ;keZ) 

(1.7) Rcu = d 2k (*gZ) 

(1.8) R r+U = d 2r+2k+2 (keZ) 

In the following, we will also consider the fully non-commutative version of (11.61 ) 
in the case of A i , namely the recursion relation 

(1.9) R„+iR„ 'R^-i =R„ + R„ 1 

for variables R„ G A, a unital algebra, and which corresponds to (11.61 ) for r = 1 but 
with different boundary conditions Ro,„ = R2, n = 1, the unit of A. 

1.2. Initial data. To fix entirely their solutions, we must supplement the systems 
above with a suitable set of initial data. 

For the Q-system, an obvious initial data set consists in fixing the values taken by 
{R a j„,R a , m +i}aei r f° r some fixed m G Z. Indeed, the system of recursion relations is 
a three-term relation in the variable k, hence fixing all R a ^ for two consecutive values 
of k allows to determine all other 7?'s. Moreover the expression for R a ^ in terms of the 
{R a ,m,Ra,m+i}aei r is the same as that for R a ,k-m in terms of x = {R a fl,R a ,i}aei r , 
by translational invariance of the equations. We call xo the fundamental initial data 
set. The most general admissible initial data sets are indexed by Motzkin paths m = 
(mi,»i2, ...,m r ) G Z r of length r— 1, with the property that —m, G {0, 1,-1}, and 
they read x m = {R a ,m a ,Ra,m a +\}aei r ( see FigH]( a )). F° r instance, the fundamental 
data set xo = x mo is coded by the null Motzkin path mo = (0,0, ...,0). The initial 
data sets x m are related via elementary moves, called mutations. We define the action 
of the (forward/backward) mutation , a G on a Motzkin path m as /i„ (m) = 
m', with m'p = mp ± 8 a p, when m' is also a Motzkin path (see FigfT] (a) for an 
illustration). Assume x m is an admissible set of initial data, then x m i differs from 
x m only by a substitution R a ,m a — > Ra,m a +2 (forward mutation /i+) or R a ,m a +\ — > 
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(a) k (b) (c) 

FIGURE 1. Initial data structure for Q-systems (a), T-systems (b) and 
the non-commutative A i Q-system ( c )- I n a h cases we indicate a sample 
forward mutation. 



Ra,m a -\ (backward mutation jx~). These substitutions are done using the <2-system 
relation for a, which involves only other terms that are common to x m and x m i. It 
is easy to see that iterated mutations allow to connect the fundamental Motzkin path 
mo to any other Motzkin path m. 

For the T-system, the obvious fundamental initial data set also consists of two 
layers with consecutive values of k, hence xo = {T a> 2j+a,Q,T a! 2j-i+a,i}jei,i U P to 
translation. The most general admissible initial data sets are indexed by "stepped 
surfaces" k = (k a j) a ei r ;j€Z with k a j € Z and the conditions that k a+l j — k a j = ±1 
and k a j + \ — k a j = ±1 (see Fig{TJ(b)). The data set associated to the stepped surface 
k is *k = {T a ,j,k a }ae/ r ;jez- The fundamental data set x = Xk corresponds to the 
surface ko with k a j = CC + j mod 2. The initial data sets xt are related via elemen- 
tary moves, still called mutations. We define the action of the (forward/backward) 
mutation j> a Ir> j € Z, on a stepped surface k as ju„ ^(k) = k', with kfp p = 
kn p ±2S a ttS p j, when k' is also a stepped surface (see Figfjlb) for an illustration). 
Assume is an admissible set of initial data, then %/ differs from x^ only by a 
substitution T a j^ aj — > T a j^ aj ±2- These substitutions are done using the T-system 
relation for a, which involves only other terms that are common to x^ and x^i. It is 
easy to see that iterated mutations allow to connect the fundamental stepped surface 
ko to any other stepped surface k. 

For the non-commutative A\ Q-system, admissible initial data have the form x& = 
(Rfc,Rjt+i), and (forward/backward) mutations lead to Xfc±i (see Figfrj(c)). 

1.3. Positivity and cluster algebra. A common feature of all these systems is that 
they share the so-called positive Laurent property, namely that their solutions, when 
expressed in terms of any of their admissible initial data, are Laurent polynomials 
with non-negative integer coefficients. 
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For the case of commuting variables this is actually part of a more general structure 
known as cluster algebras [11]. A rank n cluster algebra (without coefficients) is a 
dynamical system describing the evolution of "data" vectors x t = ((x f ),)" =1 (clusters) 
of n formal variables (cluster variables), attached to the nodes of a Cayley tree of 
degree n, with edges cyclically labeled l,2,...,n around each node. This evolution 
is coded by a companion skew-symmetric n x n matrix B t with integer entries called 
the exchange matrix, which evolves as well. The pairs (x t ,B t ) and (x u ,B u ) at two 
adjacent nodes t,u connected by an edge with label k are related via the (involutive) 
mutation fa, defined by jXk{x t ,B t ) = (x u ,B u ) with: 

J [xt)i ifi^k, 
{Xu)i ~ \ (4{n;Li(^)r d+ +nLi(^)f a]+ } otherwise. 

, . ( -(B t ) Lj ifi = koij = k; 

1 u)hJ 1 (B t )ij + sign{(B t )^)[{B t )a(B t ) k .j} + ) otherwise. 

Here we use the notation [x] + =Max(x,0). The exchange matrices B are usually rep- 
resented pictorially as quivers, namely graphs with oriented possibly multiple edges, 
with the condition that B,j =number of edges i — > j when B,j > 0. The finiteness 
restriction of rank n may be removed, while keeping the same definitions, leading to 
the notion of cluster algebras of infinite rank. The main property of cluster algebras 
is the so-called Laurent property Ifl2ll ifTTIl : every cluster variable of any node u may 
be expressed as a Laurent polynomial of the cluster variables at any other node t. The 
positivity conjecture states that this polynomial always has non-negative integer co- 
efficients jTTj . and was only proved in a few cases so far (e.g. in the so-called acyclic 
cases [3] [13] [1], or that of clusters arising from surfaces Il24l0 . 

The above Q- and T-systems have their variables belonging to specific cluster alge- 
bras in the following way |[T9lll5ll . The initial data sets turn out to be some particular 
clusters of a cluster algebra of rank 2r (A r Q-system) or of infinite rank (J-system). 
To specify the cluster algebras, we only need to give the exchange matrix Bo at the 
fundamental node of the Cayley tree, while xq = (/?io 3 ^2,Oj ••■>^r,o>^i,i> •■• I ^r,i) 
(Q-system) or x Q = {T a j^ + j mod2 ) a ei r ;jeZ (T-system) is the suitably ordered funda- 
mental set of initial data x mo or Xk defined above. For the A r Q-system, we have 

B = ^ where C is the Cartan matrix of A r , C a f, = 28 a j, — 8\ a -b\,i- F° r tne 

T-system, So is infinite. It has the entries: 

( S o)(aJ),(/U) = ( _1 ) a 1 { S j,k$\a-P\,l ~ 8 a ,p8\j-k\,l) 

with a,p el r and j,k € Z. 



1.4. Results. The aim of these notes is to present a unified approach for proving 
Laurent positivity in all three cases described above, and to extend it to more general 
cases. This is essentially based on joint work with R. Kedem [6] [7] [8] 0. The 
strategy consists of three steps. Step 1: construct quantities conserved modulo the 
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evolution equations (discrete integrals of motion), and derive a linear recursion rela- 
tion for some subset of solutions, from which all solutions can be obtained. Step 2: 
rephrase the linear recursion relation in terms of a finite continued fraction expres- 
sion for the corresponding generating function, which is manifestly positive in the 
initial data. In principle steps 1 and 2 can be repeated for each initial data set, but in 
practice, it is more efficient to only carry them out for one set of initial data (say the 
fundamental one), and then use the Step 3: express mutations as rearrangements of 
the initial continued fraction, that preserve the manifest positivity. 

Note that continued fraction language immediately translates to that of weighted 
paths on some target graphs. The Laurent positivity results may all be rephrased 
in terms of partition functions of positively weighted paths on suitable graphs. The 
mutations that implement changes of initial data (rearrangements of the continued 
fractions) may be interpreted as transformations of both the target graph and the 
weights. 

This leads to Theorem 12.41 for the A r Q-system (Section [2]), a restricted version 
of Theorem 13. II for the A r T-system (Section and Theorem 14. II for the fully non- 
commutative A i <2-system (Section RTTT) . 

This approach is actually restrictive inasmuch as T-system solutions are con- 
cerned. Indeed, it only covers the case of translationally invariant stepped surfaces 
k, with k a j + 2 = k a ,j (but arbitrary attached initial data T a jjc a ■)• The only allowed 
mutations are actually compound, i.e. infinite iterations of elementary mutations that 
map such surfaces to each-other. To include the most general case of stepped sur- 
faces, one must switch to a different point of view, by constructing for each stepped 
surface a weighted graph or network, from which the solutions can be computed ex- 
plicitly. We shall not cover this case in the present notes, and refer to JfQ for details. 

The advantage of the present approach however is that it has a natural generaliza- 
tion to fully non-commutative settings, for which such notions as non-commutative 
continued fraction rearrangements exist. We illustrate this further in Sections l4~2l and 
14.31 The general case corresponding to non-commutative weighted paths is investi- 
gated in Section l4~4l 

Acknowledgments. We would like to thank especially R. Kedem for a fruitful col- 
laboration. We also thank S. Fomin for hospitality at the Dept. of Mathematics of 
the University of Michigan, Ann Arbor, while this work was completed. We received 
partial support from the ANR Grant GranMa, the ENIGMA research training network 
MRTN-CT-2004-5652, and the ESF program MISGAM. 

2. A r Q-SYSTEM 

2.1. Warmup: the Ay case. We start with the Ay Q-system (11.11 ) for r = 1, with 

(2-1) Rn+lRn-l =R 2 „ + 1 

We have: 
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Theorem 2.1. The solution R n of the A i Q-system (12.11) , when expressed in terms of 
any admissible initial data x, = (/?,-, Ri+i), is a Laurent polynomial with non-negative 
integer coefficients. 

Proof. We wish to express the general solution R n in terms of any admissible initial 
data X, = (Ri,Ri + \). Obvious symmetries of (12.11 ) allow to restrict ourselves to R n 
for n > and xo only. Indeed, introducing the two automorphisms a,T such that 
<j(Rq) = R\ and <j(R]_) = R2, while t(Rq) = R\ and t(/?i) = Ro, we immediately 
see that o(R n ) = R n +i and z(R n ) = R\- n for all ngZ, respectively from the in- 
variances n — > n + 1 and n — > 1 — n of the relation (12- lb - If we have an expression 

= /«(xo), by applying to it & we get an expression R n = o- i '(/, 1 _ ! (x )) = g n ,i( x i) 
in terms of any x,-. Likewise, starting from R n = f n (x ), n>0, and applying T, we get 
R\_ n = t(/„(xo)) = gn(xo) which gives an expression for all /?„, n < in terms of xo. 
Moreover, if /„, n > is a Laurent polynomial with non-negative integer coefficients, 
so are all the transformed expressions under iterations of a or T. 

We now note that eq. d2.ll ) expresses that the "discrete Wronskian" 



(2.2) W n = 

Therefore, we may write 

w n+1 -w n = o 



Rn+\ Rn 
Rn Rn 1 



1 



Rn+2 + Rn Rn+l + 
Rn+1 Rn 



hence the columns of the resulting matrix must be proportional, and we write c n = 
R "+ 2 + R " = = c „_! = c, a constant independent of n. We compute: 

R n -t-] 1 Rn R] 1 R() 

(2.3) c = -2±i + h — — = — + + — 

RnRn+l Rn+\ Rq R\Rq R\ 

This gives a conserved quantity (discrete first integral of motion) for the relation d2.ll ). 
It may also be rephrased into a linear recursion relation for the /?„'s: 

(2.4) R n+1 -cR n +/?„_i =0 

Introducing the generating function Fo(f) = L,T=o f " Rn/Ro, and defining the weights 
yi =Ri/Rq, y2 = 1/(RqRi) and V3 = Rq/Ri = l/ji, we may use d2.4| ) to get: 



(2.5) 



! + 1-?(J2+B) 



1-rt + f 2 l-t(yi+y 2 +y3)+ymt 2 
1 1 



1 — i — ti — ; — t 1 



l-/(j2+y3) 1 . yi 



l-ty 3 

The latter finite continued fraction expression has manifestly positive polynomial 
coefficients of yi,yz,y3 hi its series expansion in t, and the Theorem follows from the 
fact that yi ,y2,y3 are Laurent monomials in the initial data xq. □ 
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The continued fraction expression above may be interpreted as the generating 
function for paths on the integer segment [0,3], with steps ±1, starting and ending at 

and with step weights 1 for steps i— 1 — >■ i, i = 1,2,3 and tyi for steps i — y i — 1, 

1 = 1,2,3, each path being weighted by the product of its step weights. 

The two possible forward mutations of the cluster algebra associated to (12.11 ) read: 
Hf(R 2i ,R 2i+ i) = (#2i+2,#2i+i) and fl^(R 2 i,R2i-i) = (^2i,^2i+i), where the updates 
of initial data are made by use of the relation (12- 1 b - Let us express explicitly the 
effect of the mutation /j.^ on the continued fraction expression ( 12.51) , by noting the 
two following rearrangement Lemmas, both easily proved by direct calculation: 



Lemma 2.2. 



a 1 



1 — a — b 

1 



Lemma 2.3. 



provided 



We write: 



a + ■ 



l-b 

b a' 



1 



1 — u l—c' — u 



be , ac 



a' = a + b b' = 



a+b a+b 



F (t) = =l + ty 1 F 1 (t) 



i-t- 



i-t- n 



1-^3 
1 



i-tyi-t—— l-t- 



l-ty' 3 

where we have first applied Lemma [2721 and then Lemma 1231 with u = and: 

Rf + l _R 2 



y'i = y\+y2 



}'2 



I 

>'3 



RqRi R\ 
yiyi i 



yi+y 2 R\R 2 
y\y?, R\ 



yi+y 2 R 2 

We conclude that is precisely the generating function Lnxi^^i+i/^i where 
the 7? n 's are expressed in terms of xi = jj.^(xo). Note also that F\(t) = o(Fo(t)), 
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upon defining a(t) = t. More generally, repeated application of the rearrangement 
Lemmas |2 .21 and 1231 leads to iterated mutations. 

We now have all the ingredients necessary for the generalization to the A r Q- 
system: integrals of motion, continued fractions, rearrangements. 

2.2. The A r case. The purpose of this section is to outline the proof the following: 

Theorem 2.4. The solution R a ,n of the A r Q-system is a positive Laurent polynomial 
of any admissible initial data x m , for any Motzkin path m of length r—l. 

As before, the symmetries of the system allow to restrict ourselves to R a n with n > 
and to the Motzkin paths in a fundamental domain M r modulo global translation of 
all indices by 1 . 

The A r <2-system (11.11 ) possesses r conserved quantities generalizing (12.31) . These 
are obtained by repeating the argument of the r = 1 case. 
Introducing the discrete Wronskian 

(2.6) W a ,n = det (R\, n -a-b+a+\) 

\<a,b<a v ' 

we have the relation W a , n +i^a,n-\ = W^ n + W a +\ M W a -i,n as a consequence of the 
Desnanot-Jacobi relation between the minors of any square matrix, also used by 
Dodgson for his famous condensation algorithm [ 10]. Moreover, we have W^ a =R\ >n 
and Wb,„ = 1 by definition. We conclude that W a , n = Ra,n, the solution of the A r Q- 
system. As a consequence, we have the second (A r -)boundary condition: W r +i >n = 1- 
Writing W r+ \ M+ i —W r+ \ A = 0, we deduce the existence of a linear recursion relation 

r+l 

£(-i)%/?„_ fl = o 

a=0 

with co = c r+ 1 = 1 and where c\ , C2 , . . . , c r are r conserved quantities modulo dl.ll ). In 
Ref.[6], these were expressed explicitly in terms of the fundamental initial data x m „ 
corresponding to the null Motzkin path mo, as partition functions for hard particles 
on some target graph. 

Introducing the generating function F mo (t) = Ef B /?in//?i,o> the result of [6] may 
be rephrased into the following continued fraction expression: 

(2.7) F mo (t) = 
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a+1 




a 



a-1 



(i) 



(ii) 



FIGURE 2. The two cases (i) and (ii) of forward mutations jj.^ that allow 
to reach by iteration any Motzkin path, up to global translation, from the 
null path mo. 

where the weights = y,(mo) are the following explicit Laurent monomials of the 
initial data x mn : 



Remark 2.5. The continued fraction expression (12.71 ) has an immediate path inter- 
pretation. We consider the paths on the integer segment [0, r + 2], from and to 0, 
with steps ±1,0 with respective weights I for steps i — 1 — > i, ty2ifor steps i+ 1 — > i 
(i = l,2,...,r), tyn-\ for steps i—}i((i = 2,3,...,r), andty\for the step 1 — » 0, fj2r+i 
for the step r + 2 — > r + 1. Then Fq(?) is the sum over all such paths of the product 
of their step weights. Alternatively, this was interpreted in 0] as the partition func- 
tion for weighted paths on a graph T mo with 2r + 2 vertices labeled 0, 1,2, ...,r + 
2, 2', 3', and adjacency matrix A with non-vanishing elements A^- + i = A,-+i ; = 1 
(i = 0, 1, ...,r+ 1) and Ajji = A^, = 1 (i = 2,3, ...,r). 

The expression (12.71 ) has manifestly positive polynomial coefficients of the y's in 
its series expansion in t, and the Theorem 12.41 is proved for R\ n and the initial data 
x mo . To reach all the Motzkin paths in the fundamental domain M r , we proceed by 
induction under (forward) mutation. We note that all the Motzkin paths in M r may 
be obtained from mo by iterated application of mutations }i+ : m h-> m', restricted 
only to the two following situations depicted in Fig f2] (where by convention we must 
simply omit mo and m r+ \): 

• Case (i): m„_i =m a = m a+x - 1, then m' a _ 1 =m' a -l= m' a+l - 1. 

• Case (ii):ra a _i = m a = m a+ \, then m' a _ l = m' a — 1 = m' a+l . 

In Ref.[6], it was found that the mutations may be implemented iteratively on the 
continued fraction .F mo (?) (12.71 ). to reach any F m (t) = Y*n>o tn R\,n+m { / R\,mi, m G 3VC r , 
by application of the Lemmas 12.21 and 12.31 according to the case : m i-> m' at 
hand: 

• a > 1: Apply Lemma [231 to F m (t) -» F m /(f). 

• a = 1: Apply first Lemma 12.21 to obtain F m (t) = 1 +tyi(m)G m (t), and then 
Lemma[23]to G m (f) -+F m >(t). 



(2.8) 



y2a-l — 



R<x,oRa-l,l 



RaflRa.l 
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Each continued fraction thus obtained is manifestly positive, as the two rearrange- 
ments are manifestly positive. Note that under the rearrangement Lemma 12.31 the 
structure of the continued fraction changes, and in general we obtain multiply branch- 
ing continued fractions. Remarkably however, the continued fraction F m (t) may be 
expressed entirely in terms of a collection of weights y(m) = (y\ (m), ...,y2 r +\ (m)) 
called skeleton weights. Alternatively, the continued fraction expressions allow to in- 
terpret the ratio Ri >n+m /Ri^ as the partition function for weighted paths on a rooted 
target graph r m , from and to the root, and with exactly n descents, i.e. steps taken 
towards the root. A compactified version of these graphs was found recently in H, 
and allows to put all the continued fractions F m in the form of Jacobi-type continued 
fractions, namely of the form 



Jx) 



1 



1 — Xl 



x 2 



1 -X3 



A'4 



1 -X 2r 



X2r 



and we have 

(2.9) F m (t) = l+fy 1 (m)/(y 1 (m),...,y 2r+1 (m)) 

where for a = 2,3,..., r; 

'ha-l = ty%a-\ ha = tyia if m a +l = m a 

$2a-l = t(yia-i +y2a) ha = fylayia+X if =m a + l 



3-1 " 

and hr+i = tyir+\ 



ha-i = tyia-i - &a = ^ if m a+] =m a -\ 



We see that the price to pay for having a nice Jacobi-type continued fraction structure 
for F m is that it is no longer manifestly positive, due to the subtractions in the defini- 
tion of y. However, it is a straightforward exercise to rearrange ( 12.91) into a manifestly 
positive, possibly multiply branching continued fraction, in which all the coefficients 
are Laurent monomials of the y, (m)'s with coefficient t (we call this form canonical). 
We simply apply iteratively from the bottom to the top of the Jacobi-like continued 
fraction the Lemma [2721 in the form 

be b 

(2.10) a + bA = c?H 

X—c — u 



1 



1 



or its "converse": 
Lemma 2.6. 



b - b + ^u 

- + - — = a + - 

c 1 
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Note that multiple branching will occur when applying Lemma I2T61 because the 
remaining part u of the continued fraction is now duplicated. Moreover, the iteration 
of the procedure will produce new weights, which were called "redundant weights" in 
B, that are all Laurent monomials of the skeleton weights. Note finally that we must 
apply (12.101 ) at "level zero" to F m (t) itself, with a = 0, b = \, c = y\ (m) (i.e. Lemma 
I2.2I ). The manifestly positive canonical form of F m (t) can be directly interpreted 
as the partition function of paths on a rooted target graph r m , with step weights that 
are either trivial (with value 1) or Laurent monomials of the y,'(m)'s with coefficient 
t. 



Example 2.7. The fraction corresponding to the Motzkin path (2,1,0) for r = 3 
reads: 



F 2 ,idt) = l + t W 



l-ryi + 



yi >3_ 



y 3 y± 

>5 , ye 

1-tys-t 



l-tyj 



l-t- 



yiy* , >'3.V5 
>"3 l-tyj 

yi h 



ye 

l-tys-t 



l-t- 



l-ty 7 



y4 + - 



1 - ty 7 

l-ty 3 -t 



ye 

l-tys-t 



l-ty 7 



The first expression is the Jacobi-like form. The second expression is the canonical 
form, obtained by applying Lemma |Z61 first on the third level of the continued frac- 
tion, and then on the second, and by finally applying the Lemma |Z21 It is manifestly 
positive, and introduces "redundant" weights tyiy^jyT,, ty2yAye/ iy^ys)' tyAye/y5> all 
of which are Laurent monomials of the skeleton weights {y\ ,...,y 7 ) with coefficient t. 
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Example 2.8. The fraction corresponding to the Motzkin path (0,1,0) for r = 3 
reads: 



Fn.i.o(t) = l+t- 



l-t(yi+y 2 )-t 2 - 



T4 ys 



l-tys-t 



ye 

l-ty 7 



l-t- 



y\ 



l-t- 



)'2 



l-t- 



>'3 



y4 + 



i-t- 



V42ii 

>'5 

1-tyt 



l-ty 5 -t 



ye 

l-tyj 



We have applied the Lemma \2^6\ to the third stage of the initial continued fraction, and 
then Lemma \2?2\ to its second and first stages. We get a redundant weight ty^y^jy^. 

Finally, we get the following explicit expression for the skeleton weights in terms 
of the initial data x m , for any m € M r : 

Theorem 2.9. The skeleton weights associated to the motzkin path m € M r read: 



(2.11) y 2a -i(m) 

V2a(m) 



(2.12) 



Ra,m a + \Ra-\,m a -i 
Ra,m a Ra-\,m a -i + \ 
( ^a+l,m~i i +lRa+l,m n 



-p — ifm„=m a+ i + l 

a+\,m a+l 

1 otherwise 



, Ra+l,m a + lRa-\,m 
R<x,m a + lRa,m a 
( R 2 



ifm a =m K _i - 1 
1 otherwise 



Proof. This is proved by induction on mutations. First, the formulas (12.1 1112.121 hold 
for the fundamental case m = mo (see (I2.8I >). Assume they hold for some Motzkin 
path m. As before, we only have to examine the cases (i) and (ii). As a consequence 
of Lemma l2.3l the mutation pL a : y = y(m) 14/ = y(m') acts on the skeleton weights 
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as follows: 



(2.13) 



(2.14) 



Case (i) 



Case («) 




y%a-\ +yia 

y2a+\y2a / {yia-l + J2a) 
yia+iyia-l I (j2a-l +^2a) 

yia-\ +yia 
y2a+\y2a/{yia-\ +yia) 
y2a+iy2a-l/(y2a-l +y2a) 
yta+iyia-l I 'iyia-l + J2a) 



and y'p = yp in all the other cases. Finally, it is straightforward to check that these 
relations are satisfied by the weights (12.llH2.12l ). as a consequence of the Q-system. 
The Theorem follows. □ 

Let us now consider F m (t) in its manifestly positive canonical form, with possible 
redundant weights, which are all monomials of the skeleton weights, with coefficient 
t. The series expansion in t has coefficients that are polynomials of the skeleton 
and redundant weights, hence Ri^+m, is itself a positive Laurent polynomial of the 
initial data x m . This completes the proof of Theorem 12.41 for a = 1. The case of 
a > 1 follows from a slight generalization of the Lindstrom-Gessel-Viennot theorem 
lt23l lfT7l . which allows to interpret the discrete Wronskian expression (12.61 ) for R a n = 
W a , n as tne partition function for families of strongly non-intersecting paths on the 
same target graph r m , with the same step weights (see Ref.H for details). The 
Theorem 12.41 f olio w s . 



3. A r T-SYSTEM 

In this section, we present the solution of the A r T-system when the boundary 
stepped surface is periodic. It appears to be a non-commutative generalization of the 
<2-system case, and to have the Laurent positivity property. We actually have the 
following more general result, proved in iH : 

Theorem 3.1. For any stepped surface k = (k a j)aei r ;jeZ> tne solution T a j^ of the 
A r T -system is a positive Laurent polynomial of the initial data x^. 

3.1. Periodic stepped surfaces. In this section, we restrict ourselves to initial data 
Xk, with periodic stepped surfaces k with fc«j+2 = k a j for all a,j. An example is 
provided by the fundamental stepped surface ko. These periodic stepped surfaces 
are related via compound mutations of the form jj,^ = II/ezMa /> namely infinite 
products of elementary mutations, which are implemented by the use of the T-system 
relations for a fixed a but for all j 6 Z simultaneously, in order to simultaneously 
substitute T aJMjmoi2 -> T aJMjmoi2±2 for all j G Z. 

The periodicity property allows us to work within the framework of the non- 
commutative <2-system formulation of the T-system dl.5Hl.6l ). The compound muta- 
tion indeed amounts to a substitution JH a ,k R-a,fe±2 by use of the relation dl.6l) . 
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Note that the operators R a ^ are only "mildly" non-commuting, as all the opera- 
tors of the form d^ 2a ^ 2k ^ m R a M^ m act diagonally on the basis \ j)jez, hence they all 
commute, for all a,k,m. 

3.2. Non-commutative approach. In Ref. [7], the constructions of the Q-system 
case are generalized to the case of the T-system (conserved quantities, linear recur- 
sion relations, partition functions for hard objects, and finally path formulation). In 
particular, the Desnanot-Jacobi identity may still be applied to eliminate the T a j^ for 
a > 1 in terms of the Tij^'s via the following "discrete Wronskian" formulas: 

(3.1) T a jk= det m j- a +b ,k-a-b+a+\) 

\<a,b<a v J 

We may reformulate the final result for the solution T\ of the A r T-system as a 
function of the fundamental initial data Xk as follows. The generating function 
Fo(?) = L&xi^RiitRio nas a con ti nue d fraction form analog to ( I2.7I ). Define the 
operators y a , a = 1,2, ...,2r+ 1 by: 

u\ T a t +a-i lTa-l ,t+a,0 ~\ 

y2«-il0 = -=, Tf, \t-2) 

la,t+a-lla-l,t+a,\ 

i.v 7o+l,f+a,l^a-l/+o!+l,0 1 . \ 

72a\t) = 7f ^ V- 1 ) 

*a,t+a,\l<x,t+a,0 

then we have: 

(3.2) F o (0 = 

(l-t(l-t(l-ty 3 -t(---(I- ty 2r -i - t(I - ty 2r +iy l y2 r y K ■ ■ ) ^4) "^2) ~Vl) 

where / denotes the identity operator I\j) = \j). 

Picking the coefficient of | j — k — 1 } in Fo (t ) | j + k — 1 } , this may be interpreted as 
a direct generalization of the <2-system result: T\j^/T\j + k,o is the partition function 
for paths on the same graph r mo as in the commuting case (see Remark [231 ), but now 
with "time"-dependent weights y a (f ) associated to the operators y a \t) = y a (t ) \t — 2) , 
such that the paths start at time to = j — k—l and end at time t\ = j + k—l, and that 
the "up" steps +1 take no time (i.e. do not increase the time counter) and have weight 
1, while the down steps —1 and level steps take two units of time (i.e. increase by 
2 the time counter) and receive the weights yt(t), where t is the value of the time 
counter after the step is completed. This shows explicit Laurent positivity for T\j^ as 
a function of Xk . (To make contact with [7 ], we have to introduce conjugated weights 
Y a = d~^~ l y a d^ +1 for a > 2 and Yi = dyid -1 , with the effect of distributing the 
time more evenly on the steps, namely steps ±1 take one unit of time, while steps 
take two.) 

Finally, let us show that, like in the commutative Q-system case, we may imple- 
ment mutations by non-commutative versions of the rearrangement Lemmas 12.21 and 
12.31 We have the following two non-commutative rearrangement of non-commutative 
fractions, which preserve manifest positivity. (In the following, and for later gener- 
alizations, we assume the continued fraction coefficients are elements of an algebra 
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A with unit 1; moreover, any expression of the form (1 — a) -1 may be formally ex- 
panded as L„>oa"-) 

Lemma 3.2. Let a,b € A, we have: 

l + (l-a-b)- 1 a= (l-tl-b)-^) -1 

Lemma 3.3. For a, b, c, u € A, a + b invertible, we have: 

a+(l-(l-u)- 1 cy l b=(l-(l-c'-u)- 1 b , y l a' 

provided 

a' = a + b b' = cb(a + b) _1 c , = ca(a + b)" 1 

We note that the periodic stepped surfaces k may be indexed by Motzkin paths 
of length r — 1. Indeed, the map <p : k i-» m defined by m a =Min(k a fi,k a i) is a 
bijection between the set of periodic stepped surfaces and that of Motzkin paths of 
length r— 1. By a slight abuse of notation, we will denote by F m (?) the function 
Fk(t) with m = <p(k). Starting from the operator continued fraction F mo (t) (13 -2b . 
we may now apply iteratively the above rearrangement Lemmas 13-21 and 1331 exactly 
in the same way as we applied Lemmas 12.21 and 1231 for the commutative Q-system. 
The results are operator continued fractions F m (t), with operator coefficients, whose 
matrix elements in the basis \j)jeZ are Laurent monomials of the initial data x^ with 
coefficient t. As the rearrangements preserve positivity manifestly, we deduce that 
the Theorem 13.11 holds for a = 1 and all periodic stepped surfaces k. Alternatively, 
the rearranged continued fractions F m (t) generate time-dependent weighted paths on 
the same target graphs r m as in the commuting Q-system case. 

For a > 1, we interpret the discrete Wronskian formulas (13.11) via the same gener- 
alization of the Lindstrom-Gessel-Viennot theorem, as computing the partition func- 
tion for families of strongly non-intersecting paths on r m , now with time-dependent 
weights. The positivity and therefore Theorem 13. II follows, for the case where the 
stepped surface k is 2-periodic in the j direction. 

4. Non-commutative systems 

We have already seen how to interpret the T-system as a non-commutative Q- 
system, in terms of variables obeying many special commutation relations. In this 
section we investigate integrable relations involving fully non-commutative variables. 

4.1. The non-commutative A i <2-system. We have the following: 

Theorem 4.1. The solutions of the non-commutative A\ Q-system d 1 -91) in terms of 
any admissible initial data x, = (R,-, R;+i ) are Laurent polynomials with non-negative 
integer coefficients. 

To prove the Theorem, we must somehow repeat the steps of Section 12.11 As in 
the commuting case, the symmetries of the equation allow us to restrict to the case 
n > and to the initial data Xq. We still define anti-automorphisms a, z that leave 1 
invariant, and such that <p{ab) = q>{b)(p(a), <p = <t,t, and cr(Ro) = Ri, cr(Ri) = R2, 



16 



PHILIPPE DI FRANCESCO 



while t(R ) = Ri and t(Ri) = Ro- We deduce that a(R„) = R„ + i and t(R„) = Ri „ 
for all Assuming we have a positive Laurent polynomial expression for R„(xo) 

for all n > 0, we may obtain a positive Laurent expression for n < as well by 
applying T. Finally from a positive Laurent polynomial expression for R„(x ) for all 
n S Z, we may obtain one for R„(x,) for any i G Z, by applying a 1 . 
We now show the existence of conserved quantities: 

Lemma 4.2. The equation d 1 .9b has the two following conserved quantities: 

(4.1) C = R^R^+iR- 1 

(4.2) K = (R^+CR^OR-^R-^+iC + R^i) 

Proof. Introducing C„ = R^RnRn+iR" 1 , and noting that R„ commutes with the 
r.h.s. of (11.91 ), we have R,^ + 1 = R M+ iC M _iR„_i = R M+ iC M R„_i, hence C„ = C M _i = 
C. Note that the conservation of C„ implies a quasi-commutation relation: 

(4.3) R M R,j+i = R„+iCR M 
and that the A\ Q-system may be rewritten as: 

(4.4) R„+iCR„ i =R" + 1 

Introducing K n = (R„ +1 +CR„_i)R" 1 and L„ = R-^R^+jC + R„_i), we find that 
R„K„R„ = R„L„R„ by use of the quasi-commutation (I4.3I ). hence K„ = L„. More- 
over, computing 

R«+iK„R„ = R,^ +1 +R,^ + 1 = R„+iL„ + iR„ 
we deduce that K„ = L, !+ i = K, !+ i = K, and the Lemma follows. □ 
The conserved quantities may be rewritten as: 
C = R„ H 1 1 R„R, 1+ iR„ 1 =y3yi 
K = R„ +1 R,7 1 +R^ 1 R; 1 +R^ 1 R„ = y 1 +y 2 + y 3 
where we have defined: 

(4.5) y^RiRo 1 y 2 =Rr 1 R 1 y 3 = Rr 1 Ro 

The Lemma implies the existence of a linear recursion relation: 

R/j+i — KR„ + CR„_i = 
For t a central element in A, introducing the formal generating function Fo(t) = 
Ln>o f "R«Ro ' we easi ly compute: 

F o (0 = (i-(yi+y2 + y3)? + y 3 yi? 2 ) _1 (i-Ky2 + y3)) 
= ^i- f (i- f (i- fy3 )- 1 y2 ) _1 yi 

This may be expanded as a series in t, whose coefficients are polynomials of yi ,y2,y3 
with non-negative integer coefficients. The theorem |4~T1 follows from the definition 
of Fo(t) and the explicit form of the y's (14.5 b - 
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Finally, let us show that, like in the Q and T-system cases, we may implement 
mutations by use of the non-commutative rearrangement Lemmas 13.21 and 13.31 We 
write: 

F (0 = (l-f(l-f(l-?y 3 )- 1 y2) _1 yi) _1 =l + ^F 1 (0yi 

Fi(0 = (i-<yi-*(i-<y3)- 1 yi)" 1 = (i-t(i-(i-y , 3rV 2 )"Vi)~ 1 

where we have first applied Lemma [3721 and then Lemma [331 with u = and: 

y'i = yi+y 2 = (Ri+Rr 1 )R 1=R 2Rr 1 
y' 2 = y 3 y2RiR 2 1 =Rr 1 C R^R.Rr 1 
= y 3 yiRiR 2 1 = CRiR 2 1 =R 2 1 R 1 

where we have used the Q-system relation dl.9l > (first line) and the quasi-commutation 
relation (14.31) (second and third). 

We conclude that Fj (t) is precisely the generating function £„> f'^+iRj" 1 where 
the R„'s are now expressed in terms of Xi = /ij + (xo). More generally, repeated appli- 
cation of the non-commutative rearrangement Lemmas 13.21 and 13.31 leads to iterated 
mutations. 

4.2. The affine non-commutative rank two cases. The rank two affine cluster al- 
gebras have a fundamental skew-symmetrizable exchange matrix So = 

with (b,c) = (2,2), (1,4) or (4, 1). The A\ <2-system of previous section reduces to 
the case b = c = 2 when C = 1. The two other cases also have non-commutative 
counterparts, namely the system 

R2nR 2 „ 1 _iR2«-2R2n-l = 1+R 2 n-1 

(4.6) R 2 „+iR 2n 1 R 2 „-iR2, 1 = l + (R 2 n) 4 

and that obtained by interchanging even and odd indices. 
We have the positivity Theorem: 

Theorem 4.3. For all n € Z, the solutions R„ of (14.61 ) are positive Laurent polyno- 
mials of any initial data x,- = (R,,R 1+ i). 

Repeating the approach of previous section, we note that the symmetries of the 
system allow to restrict to the case of n > and to the initial data xo and Xi . Next, we 
identify two conserved quantities: 

C = R n+1 R„R„ + iR„ 1 

K = (R2n+2 + CR 2 „_ 2 )R 2f? 1 

We finally arrive at the following continued fraction expression for the generating 
function Fo(f) = L«>o f "R2«Ro 1 in terms of xq: 

(4.7) F(t) = (l-Kf + CfVtl-KK-RaRo- 1 )) = (l -; yi -t 2 (l - m )- x y 2 )~ l 



U o)' 
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where 

yi = RaR^Ml+R^R^RiRo 1 

y 2 = (K-yi)yi-C 

= (l + Rr^Ro^l + Rr^Ro^RiRo^+Rr^Rr^o'RiRo 1 
y 3 = K-^CR^RoMl + Rr^Ro' + Rr'Ro 

This proves positivity of all R 2 „ in terms of x . The case of X\ is obtained by muta- 
tion/rearrangement of Fo(f). We write: 

r l (F Q (t)-l)y^ = (l-Kf + CfVCl-tCyr 1 ) 
= (l- ty [- t i(l-tj 3 )-y 2 )- 1 

where 

yi = K-y^=R4R 2 1 =R 2 1 (l+Ri)R2 1 Rr 1 +R^Rr 1 
y 2 = y^(K-y^)-C 

= r 2 - 1 (i+Ri)R2 2 ( 1 + r i) r 2 Rr+Rr 1 

y' 3 = Cyr 1 = CR ( ,R 2 1 =R 2 1 (l + Ri)R 2 1 

We get a series in t with coefficients that are polynomial in y[ , y 2 , y' 3 , and as each of 
them is itself a positive Laurent polynomial of the initial data we deduce positivity 
for all R 2 „ in terms of Xi. The case of the variables R 2 «+i = R 2n+2 CR 2 „ — 1 is easy, 
by directly checking that R 2;1+2 CR 2 „ contains the term 1. This completes the proof 
of Theorem l4.3l 

4.3. The non-commutative A 2( t systems. We now consider the systems 

(4-8) U 2w+2( t + iU 2 „ = U2n+lU2n+2fc + 1 

(4.9) U2n+ll2n+2&+2 = U 2 , 1+2j t + iU 2 „ +2 + 1 

with initial data x, = (u,,u, + i , . . . ,u, +2j t) for all i € Z. The commutative version of 
(I4.8H4.9I) corresponds to mutations within the affme A 2 £ cluster algebra. 
We have the following positivity result (||9"T0: 

Theorem 4.4. The general solution u„ to the system A4.8M.9i is a positive Laurent 
polynomial of any initial data Xi,for all n,i £ Z. 

This is proved in the usual way. We may restrict ourselves to n > and to the 
initial data xo. We find the conserved quantity 

K = (u 2 „ +2 + u 2m _ 2 )u 2 ,/ =u 2n 1 +1 (u 2 „ + 3 + u 2 „_i) 

k 

= U U- k l + U 2i U 1 + £ U 2; ^ ! (U 2 / + U 2; ^ 2 ) 

7=1 
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This allows to derive continued fraction expressions for F,(?) = L/!>o f " u 2i+2fai and 
=Lr=o f " u 2i+i+2/tn, i = 0, l,...,k- 1: 

F,-(0 = (l-(y 2 i + Z2i)t + t 2 )- l (l-Z2it)U2i 

= (i - ynt - (i - Z2ity l w 2i t 2 y [ u 2i 

Gi(t) = U 2;+ l(l-Z2,- + lO(l-(y2 I -+l+Z2 ( -+l)? + ? 2 ) _1 

= u 2i -+i (l-/y 2i+ i -f 2 w 2 ,- + i(l-z 2! - + i0~ 1 ) 1 

where 

f2i = U2,+2,tU 2 / yu+l = U 2 ,-+iU2,-+l+2/t 

Z 2 ; = K — y 2 ; = U2 ; --2/tU 2! - 1 Z2i+1 = K — y 2; -+i = U 2( . I | _ 1 U2, + l-2/t 

W 2i ' = Z 2 ;y 2i ' - 1 W 2 i + l = y2*+l z 2i+l - 1 

We have manifest positivity in the y,z, w's. Moreover, the y's are defined recursively 
by yo = u 2 *Uo 1 and y 2; - = y 2 j-\ +u 2 j_ ] u 1 j, y 2j+ i = yij+f^j^vqj, hence they form 
subsums of K, and so do the z's. Finally all y's contain the term U2&UQ , while all z's 
contain UqU^V hence both y ; z, and z ; y ; - contain 1: we conclude that the y,z,w's are 
all positive Laurent polynomials of xo. This completes the proof of Theorem l4.4l 

4.4. Non-commutative paths: a general study using quasi-determinants. In view 
of the accumulated evidence for path models and continued fractions to play an im- 
portant role in the discrete non-commutative integrable systems, we now address the 
general case of paths on the graphs r m of the A r <2-system solution, but with fully 
non-commutative weights. 

For any m Motzkin path of length r — 1 , we define the family of non-commutative 
continued fractions F m (f) as follows. Let (yi(m),...,y2 r+ i(m)) G A 2r+l be a given 
collection of non-commutative "skeleton" weights. We define the non-commutative 
Jacobi-type fraction J(x) for a family of non-commutative weights x = (xi , X2 r+ 1 ) € 
yi 2r+1 by: 
(4.10) 

J(X) = (l-X 1 - (I-X3 - (...(1-X 2r _l - (l-X2 r+ l)~ 1 X2,.)~ 1 ...)~ 1 X4)~ 1 X 2 ) 

We set 

(4.11) F m (0 = l+?J(yi(m),...,y2 r+ i(m)) yi (m) 

where for a = 2,3, ...,r: 

yia-\ = ty 2a -\ ha = ty 2a if «?a+l = m a 

ha-l =t{y 2a -\+y2a) f2a = t 2 y2a+iy2a if m a+ l = m a + 1 

ha-l = ty 2a -\ - y2a + iy 2 « = y2a+iy2a if m a +l = m a - I 

and y 2 ,-+i = y 2r +i 

Like in the commutative case, we may easily bring the continued fraction F m (f) to 
a manifestly positive form where all non-trivial coefficients are Laurent monomials 
of the y's with coefficient t (we still call this the canonical form), and for which the 
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coefficients of the formal series expansion in t are therefore Laurent polynomials of 
the y's with non-negative integer coefficients. This is readily performed by applying 
iteratively from the bottom to the top of the Jacobi-type continued fraction the non- 
commutative rearrangement Lemma 13^21 in the form 

(4.12) a+b + (l-c-u) -1 cb = a+(l-(l-u) _1 c) _1 b 

or its "converse": 

Lemma 4.5. For all a , b e . u 6 A, c invertible, we have: 

a-c" 1 b + (l-c-u)" 1 c" 1 b = a+(l-c-u)" 1 (b + uc" 1 b) 

The net effect is to "undo" the y's, at the expense of creating branchings and new 
"redundant" weights that are Laurent monomials of the skeleton weights, all with 
coefficient t, and therefore all manifestly positive. 

The manifestly positive canonical form of F m (t) can now be directly interpreted as 
the partition function of non-commutative weighted paths on the same rooted target 
graph r m as for the commuting case. 

The continued fractions F m (?) are related via non-commutative mutations as fol- 
lows. Restricting like in the commutative case to the cases (i) and (ii) for m and 
m' = /x+(m), we see that the forward mutation /x+ is implemented on F m (t) written 
in the Jacobi form (14.111 ) via the following sequences of rearrangements: 

• Case (i):(a) Undo y2«-i,y2a by applying Lemma l3T2~l (b) Apply Lemma [331 
to the piece containing y2«-i,y2a and y2a+i in factor (c) Go back to the 
Jacobi form by using (14.121 ) in reverse 

• Case (ii):(a) Apply Lemma|33]to the pie ce con taining y 2 a-i,y2a,y2«+i,y2a+2 
(b) Go back to the Jacobi form by using (14.121 ) and Lemma 1431 in reverse 

This produces F m i(t), in Jacobi form as well. Note that when a = 1 we have to apply 
(14.121 ) at the level zero of the continued fraction, with a = 0, b = 1 and c = yi (m). 
So the transformation of the continued fraction F m (?) reads: 

a = l: F m (0 = l + fFtf(0yi(m) 
o>l: F m (0 = F m ,(t) 



It is easy to follow the transformations of skeleton weights in the process. The 
mutation pL a : y = y(m) (->• y' = y(m') acts on the skeleton weights as follows: 

'y'la-l = 



(4.13) 



(4.14) 



Case (i) 



Case (ii) 



y'la 

y2a+l 

y'za-i 
y'ia 

y2ce+l 
,y2a+2 



y2a-l +y2a 
y2a+iy2a(y2a-l +y2a)~ 1 
y2a+iy2a-l(y2a-l +y2a)~ 1 

y2a-l +y2a 
y2a+iy2a(y2a-l +y2a)~ 1 
y2a+iy2a-l(y2a-l +y2a)~ J 
y2a+2y2a-l(y2a-l +y2a)~ 1 
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while all other skeleton weights are left unchanged. 

We now wish to relate the skeleton weights y a (m) to the coefficients of the series 
expansion for F m (f). In particular, we set F mo (f) = £ n>0 £ n R n RQ . 

As shown in @ this can be done by the use of the theory of quasi-determinants 

ED ED. 

For any square kxk matrix A = (a,j) with entries in A, and any p,q G {1,2, ...,k}, 
we define the (p,q) -quasi-determinant |A| M as: 

\Mp-q = a p,q~ XL a Pj(l^lij) & h<l 

whenever the r.h.s. is well-defined. The quasi-determinant has many interesting prop- 
erties (see |fl5ll for a detailed study). In the commutative case, it reduces to a ratio 
of determinants, namely \A\ p ^ q = \A\/\A p ' q \, where A p - q stands for the matrix A with 
row p and column q erased. We now define the quasi-Wronskians of the sequence 
(R n )nez to be: 



(4.15) 



Rn+a-1 
R«+a-2 

R;i 



Ri+ce-2 
R/j+a-3 



R„ 

Rf7 1 



R 



n-1 



R 



n-a+l 



1.1 



with the convention Ao,« = 1. These quasi-Wronskians satisfy a generalized Desnanot- 
Jacobi relation (also called the non-commutative Hirota equation ||9ll): 

( 4 -!6) A« + i,„ = A«,„ +1 - A aM (A a ) 1 _ l -A^_ l n )A a ^ 

The A r boundary condition is replaced by 

A , M = 1 A r+2 ,„ = (n G Z) 

expressing simply that the continued fractions are actually finite (y2r+2(mo) = 0), or 
in other words that there exists a vanishing left linear combination of the R„'s with 
r + 2 constant coefficients (the integrals of motion), causing the quasi-determinant to 
vanish (this linear recursion relation for R„, n > is then used to extend the defini- 
tion of R„ to n < as well, as the unique solution to the linear recursion relation). 
Introducing new variables 



ce-l 



(4.17) 



D a ,n — Aa+i,/ I A a n C a ,n — \\ Ace-i,«-iA a _ ! - n _ i _j 



i=0 



we may recast the above into the relations: 

(4-18) D an+ i = C a+ \ n+ iC a ^ n D an C a -\ n C a ^ n+l 

(4.19) C a ,„+iC a _ ljJ = D an + C a ^C a ] i ln 

(4.20) C a+u C a 

Using these relations, we get the non-commutative version of d2.llH2.12l : 
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Theorem 4.6. The skeleton weights for F m (t) are expressed in terms of the R„ 's as: 
(4.21) y 2 «-i(m) = C a , ma+ iC a i lma i+1 

' C a+hnia+l+1 C a l +l +1 if ffl B =m«+i + 1 



y2a(m) 



1 otherwise 



(4.22) x ^xf^ 1 ^' ^'"otherw'ise 1 " 1 } 

Example 4.7. T/jg weights for the fundamental Motzkin path mo = (0, ...,0) read: 



y 2a _i(m) =C aA C a _ 1 j y2a(m)=Z) ai 



Those for the "ascending" Motzkin path mi = (0, 1,2, ...,r — 1) are: 

y2a-l(mi) = C ce;Ce C a i lce _ 1 = A 0!!0( A~ a _ 1 

y2o(mi) = D aa = A a+ i, a A a ^ a 

This is nothing but the non-commutative version 11161 of the Stieltjes formula 11271 
expressing the coefficients of the ordinary continued fraction expansion in terms of 
those of the series expansion. In a sense, our construction is a double generalization 
of that formula. 

Note that, due to the explicit value of yi(m) = C\ m +\ = Ri jOTl+ iItf , we have 
Fm(0 = L«>o f "^i,«+m 1 Rr,i, 1 ■ We conclude that for all Motzkin paths m of length 
r— 1, the quantity R„ +W1 R~ is the partition function for non-commutative weighted 
paths on r m . All the weights are monic Laurent monomials of the skeleton weights 
of Theorem I4.6L hence the result is a positive Laurent polynomial of the variables 
C,D occurring in Theorem 14. 61 Note that these are themselves Laurent monomials of 
certain quasi-Wronskians (14. 15b - 

Except for the point mi of Example |4.7 1 above, for which we have a clear Laurent 
positive phenomenon in the variables (A ai(X , ^a,a-iYo!^v tne classification of rele- 
vant initial data in terms of the variables A an for the equation ( 14. 161 ) is not clear 
at this stage. All we have at our disposal is Laurent positivity in terms of the non- 
commutative variables C,D. A better understanding of what the "good" dynamical 
variables are should lead to a notion of non-commutative cluster algebra, via the mu- 
tations thereof. 

Apart from the cases already treated in these notes, all of which may be viewed as 
particular cases of this section, an important application concerns the so-called quan- 
tum cluster algebras 0. These are non-commutative versions of cluster algebras, in 
which the variables obey ^-commutation relations. A consequence of our construc- 
tion is that the positive Laurent phenomenon holds for all R„ = Ri 5 „, n > for the 
corresponding A r quantum <2-system Q : 

q<x(r+l a ^'R a n+l 'R a n _ l = (R a ^) 2 J t -'R a+ i n R a _i n 

with the boundary condition Ro,„ = R r +i,„ = 1, and where the variables within the 
same cluster (or initial data set) have the ^-commutation relations for 1 < a < j8 < 
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r+h 

Theorem 4.8. As a function of the initial data x m = {R a ,m a > Ra,m„+i }aei r > the quan- 
tity Ri n+mjRj - ^ is a Laurent polynomial with coefficients in Z + [q,q~ 1 ]. 

Proof. We simply have to compute the quasi-Wronskians: 

A _ -^(a-l)(2r+2-a) R R -l 
^a,n—q K a," K a-l,«-l 

and the corresponding values of the C,D variables: 

Ca,n = q 2 ^-a.nR-a.n-l D an = Ro+l,(t^a,n^a,n-l^a-l,n-l 

which are all Laurent monomials of the initial data, with possibly some integer power 
of q in factor. The Theorem follows immediately from the expressions of Theorem 
I4~6l ' □ 
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